B-mode polarization of the CMB from the second-order photon quadrupole 
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ly-^ ■ We study a new contribution to the polarization of the Cosmic Microwave Background induced 

at the epoch of recombination by the second-order quadrupole moment of the photon distribution. 
At second order in perturbation theory the quadrupole moment is not suppressed by the inverse of 
the optical depth in the tight coupling limit, as it happens at first order in perturbation theory. We 
^ ' concentrate on the B-mode CMB polarization and find that such a novel contribution constitutes 

\^ ' a contamination in the detection of the primordial tensor modes if the tensor to scalar ratio r is 

00 . smaller than a few xlO~^. The magnitude of the effect is larger than the B-mode due to secondary 

CO ' vector/tensor perturbations and the analogous effect generated during the reionization epoch, while 

CO ' it is smaller than the contamination produced by the conversion of polarization of type E into type 

, B, by weak gravitational lensing. However the lensing signal can be cleaned, making the secondary 

modes discussed here the actual contamination limiting the detection of small amplitude primordial 
' gravitational waves if r is below ~ 10" 

p ^. PACS numbers: 98.70.Vc,98.80.Es,98.80.Cq DFPD 07/A/05 

o 

a ■ I. INTRODUCTION 
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The hunting for the polarization of the Cosmic Microwave Background (CMB) is one of the major targets of 
present and future planned experiments since it offers a unique opportunity to gather information about both the 
' early (inflationary) universe and more recent epochs. The polarization is generated by the anisotropic scattering of 
CMB photons off free electrons. The decomposition into E and B-modes [l[ [1] is extremely useful when dealing with 
the nature of the cosmological perturbations that are responsible for the temperature and polarization anisotropics of 
the CMB: the B-modes are excited only by vector and tensor perturbations, thus allowing eventually to disentangle 
their effect from that of the scalar perturbations. One of the fundamental predictions of inflationary models is the 
generation of a stochastic background of gravitational waves (tensor perturbations of the metric) which makes the 
search for the B-mode polarization crucially linked to a probe of the inflationary scenario. The amplitude of this 
background is determined by the energy scale of inflation, which can widely vary among different inflationary models. 
From this point of view, future satellite missions, such as Planck, will have enough sensitivity to either detect or 
constrain the B-mode CMB polarization predicted by the simplest inflationary models. A lot of experimental and 
theoretical issues are currently under investigations in order to assess the feasibility of a detection of the primordial 
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gravity waves via the B-mode polarization. The main problematic aspects come from foregrounds 's'], instrumental 
noise Q, and the gravitational lensing on the CMB by the matter distribution. These effects might actually mask 
the signal due to primordial tensor modes. As far as the lensing "contamination" is concerned, this is a non-linear 
effect that implies the transformation of E-modc into B-mode polarization ^ and it has been pointed out that the 
inflationary gravitational-wave background can only be detected by CMB polarization measurements if the tensor 
to scalar ratio r > 10^^, which corresponds to an energy scale of inflation larger than 3 x 10^^ GeV [H, 0, B 0]- 
However, a better technique to clean polarization maps from the lensing effect has been proposed, which would allow 
tensor-to-scalar ratios as low as 10~^, or even smaller, to be probed p^Q. .11]. 

In this paper we point out the existence of a new contribution to the B-mode polarization. It is due to the 
observation made in Ref. [T^ that at second-order in perturbation theory the quadrupole moment of the photons is 
not suppressed by the inverse of the optical depth in the tight coupling limit, as it happens at first order in perturbation 
theory. The anisotropic Thomson scattering of the photon quadrupole leads to a polarization of the CMB. The study 
of the temperature anisotropies at second-order have been analyzed in several works [31 and recently Refs. [l^, [i3| 
studied the evolution of the temperature anisotropies in the tight coupling regime through the full derivation of 
the Boltzmann equations for the photon, baryon and cold dark matter components. In Ref. [T^l it has been shown 
how B-mode polarization is produced during the last scattering and reionization epochs by the non-linear evolution 
of cosmological perturbations which generates unavoidably vector and tensor modes, starting from the primordial 
(scalar) density fluctuations. 

In our case it is the non-suppression of the second-order quadrupole in the tight coupling regime that makes the 
contribution to the B-mode polarization considered here of a different origin than that of Ref. (isj . Our findings 
show that the corresponding level of the B-mode is larger and also that it overcomes by two orders of magnitude 
the analogous signal generated during the reionization epoch p^ . It turns out that such a novel contribution to the 
CMB B-mode constitutes a background to the detection of the primordial gravitational waves effect if the tensor to 
scalar ratio is smaller than a few xlO~^. The magnitude of the effect is smaller than the contamination produced by 
the conversion of polarization of type E into type B, by weak gravitational lensing. However the lensing signal can 
be cleaned, making the secondary modes discussed here the actual background contaminating the detection of small 
amplitude primordial gravitational waves. 

These secondary effects always exist and their amplitude has a one-to-one relation with the level of density pertur- 
bations, which is severely constrained by both CMB anisotropy measurements and Large-Scale Structure observations. 
Therefore, their properties are largely inflation model- independent, contrary to primary tensor modes whose ampli- 
tude is not only model-dependent, but is well-known to be suppressed in some cases, like e.g. in the so-called curvaton 
model for the generation of curvature perturbations (l7| . 

The plan of the paper is as follows. In Sec.|TT]we recall the origin of the quadrupole moment in the tight coupling limit 
at second-order in perturbation theory. In Sec. lIIII the contribution to the power spectrum of the B-mode polarization 
from the second-order quadrupole is computed and finally in Sec. IIVI we compare the results to the contribution to 
the B-mode from the inflationary tensor modes and the gravitational lensing and with other secondary effects as well. 



II. SECOND-ORDER QUADRUPOLE MOMENT OF PHOTONS AND THE TIGHT COUPLING LIMIT 

In this section we recall the starting point of our paper: at second (and higher) order in the perturbations the 
quadrupole moment of the photon distribution cannot be neglected even in the limit of tight coupling. At linear order 
the quadrupole and higher order moments of the photons, in the tight coupling limit, are negligible with respect to 
the first two moments (the energy density and velocity) because they turn out to be suppressed by increasing powers 
of the inverse optical depth t^^. However in the non-linear regime there is a source for the quadrupole made up 
of linear-velocity squared terms which do not vanish in the tight coupling limit. Such a result has been obtained in 



details in Ref. 1^, and in the following we limit ourselves to give some definitions and summarize the main steps 
necessary to achieve it. 

The quadrupole moment of the photon distribution is defined as 

n;^ = / § (nV -i5-)A, (1) 

where A — A^^^ + A'--^^ /2 are the photon temperature anisotropies split into a first and a second-order part. They 
are given by 
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with rj the conformal time, and the zeroth-order (Bose-Einstein) value around which the distribution function 
of photons is expanded, / — /^^^ + Sf — /(°^ + /^^^ + /^^V^- Eq. ^ represents the photon fractional energy 
perturbation (in a given direction), which is the integral of the photon distribution function perturbation Sf over the 
photon momentum magnitude p {p^ — pn'). The angular dependence of the photon anisotropics A can be expanded 
as 

i 1-^ 

A(x,n)=^ ^ A,M{-^YJ-^YUn), (3) 



with 



A.™ = / dnAY;„,in) . (4) 

The Boltzmann equations up to second-order for the photon temperature anisotropics A ( and the hierarchy equations 
for Aim) have been obtained in Ref. together with the second-order Boltzmann equations for the baryons and 
cold dark matter. One can integrate the R.H.S. of the Boltzmann equation for A^^^ over dflnin^n^ — (5*-'73)/47r and 
set it to be vanishing in the limit of tight coupling. Then the R.H.S. becomes 



(5) 



where is the photon-baryon velocity, S\, the baryon density contrast, \1/ the gravitational potential, and the differential 
optical depth r = fieaTCi sets the collision rate in conformal time, with n,, the mean free electron density and ctt 
the Thomson cross section.^ Therefore, in the limit of tight coupling, when the interaction rate is very high, the 
second-order quadrupole moment is given by 

n;^" ^ ^ (^vv - i<5*^v) , (7) 

as it it follows by requiring that Eq. ([5]) vanishes (the term multiplying (Sh + ^) goes to zero in the tight coupling limit 
since it just comes from the first-order collision term). At linear order one would simply get the term {—9flll^/10) 
implying that, in the limit of a high scattering rate f, Ut^ goes to zero. However at second-order the quadrupole is 
not suppressed in the tight coupling limit because it turns out to be sourced by the linear-velocity squared. 

Indeed the result of Eq. ((T]) is not surprising and had to be expected. First, the limit of tight coupling corresponds 
to treat the photons in the fluid approximation, and, at second-order, the energy momentum tensor of a fluid is 
exactly characterized by an anisotropic stress given by Eq. ([7]) (see, for example, Ref. [2^) Moreover it is also known 
that, besides the usual dipole contribution, beyond linear order a quadratic Doppler effect gives rise to a quadrupole 
anisotropy (see Refs. [H, tUl [13] ) . 



III. CMB POLARIZATION 



Since the polarization of the CMB is generated by the scattering of the quadrupole anisotropics by free electrons, 
armed with Eq. ([7]), we will explore the effects of such a contribution to the quadrupole on the CMB polarization.^ 
Notice that the resulting polarization power spectrum has a different origin from that computed in Ref. [l5l | . Ref. 
shows how second-order vector and tensor perturbations can produce CMB polarization through the decoupling phase. 
Once the vector and the tensor perturbations are obtained at second-order as a non-linear combination of first-order 



^ The perturbed line element around a spatially flat FRW background is taken in the so-called Poisson gauge [Tsl , Il9ll , which eliminates 
one scalar degree of freedom from the gm component and one scalar and two vector degrees of freedom from gij , thus reducing at linear 
order to the Newtonian gauge 

ds'^ = a^{r)) {- (1 + 21-) dr?2 - 2Vidr]dx^ + [(1 - 2$) Sij + 2Hij] dx^dx^} . (6) 

Each perturbation can be decomposed into a linear and a second-order term, e.g, ^ = -(- VI/(2)y2. Vi contains scalar and vector 

(divergence-free) perturbation modes, while Hij is a tensor mode (i.e. traceless and divergence-free, d'Hij = H'^ = 
^ The effects of the quadrupole moment jTj on the CMB temperature anisotropies have already been studied in Ref. 
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scalar perturbations, the induced power spectra of the E- and B-mode polarization are computed using the standard 
techniques developed at linear order [2^: a quadrupole is generated by the vector and tensor perturbations only 
beyond the leading order in the tight coupling expansion (otherwise it is suppressed) and to evaluate it one has to 
carefully account for the epoch around decoupling. The result is thus proportional to At;*, the width of the visibility 
function (i.e, the thickness of the last scattering surface). The case of Eq. ([7]) is completely different in that it is 
a quadrupole that is indeed generated in the tight coupling limit, and the corresponding power spectra of CMB 
polarization will be computed following different approximations as well. 



A. Polarization angular power spectra 

In order to study the generation of polarization anisotropies we will employ the formalism of the total angular 
momentum method [23 . which has the advantage of putting scalar, vector and tensor perturbations on an equal 
footing. This is the most convenient approach for our calculations, and we are going to extensively use the results of 
Ref. [24|. The temperature and polarization fluctuations are expanded in normal modes that take into account the 
dependence on both the angular direction of photon propagation n and the spatial position x, j,G'™(x, n) 



e(x,n,77) = 



_fk_ 



2 



e m=-2 
2 



(Q±*C/)(x,„,,y) ^ [A^Y^Y: (<'±.i?^^)±2Gr, 

^ ' 1 m=-2 



(8) 



with spin s = describing the temperature fluctuation and s — ±2 describing the polarization tensor and m = 
0, ±1,±2 denoting scalar, vector and tensor perturbations, respectively, i?^™^ and _B^'"'' are the angular moments of 
the electric and magnetic polarization components and 



,Gr(x,n) = H) 



2e + r 



,rr(n)]exp(*k.x) 



(9) 



Before keeping on recalling the main results of the total angular momentum method, a comment is in order here. 
The Boltzmann equations for the temperature and polarization anisotropies we are going to write are derived within 
linear perturbation theory. Then, given the second-order quadrupole ([7]), our approach will be that of evolving 
it linearly into Eqs. (|17p . In fact this is the correct procedure if we wish to isolate its contribution to the CMB 
polarization at second-order, distinguishing it from other unavoidable corrections to the equations (jl7p coming in the 
form of first-order squared terms. In a similar way, the variable A^^^ (x, n, rj) will be taken to be 40(x, n, 77), a relation 
holding at the linear level (see Refs. [T^. [26j). 

The Boltzmann equation describing the time evolution of the radiation distribution under gravitation and scattering 
processes can be written as a set of evolution equations for the angular moments of the temperature, O^™'' (for £ > m). 



and both polarization types, E^"^^ and B^"'' (for £ > 2 and m > 0) 



(m) 



e 



(m) 



{2£-l) 



(2^ + 3)^+1 



fe 



,(m) 



(m) 



B 



(m) 



(2£- 1 



■E, 



(m) 



2m 



■B 



(m) 



£{£ + ! 
2m 



-B 



(m) 



{2£ + iy 



-E, 



(m) 



2^£+l n(m) 
(2^+3)^+1 



tB 



(m) 



{2£-iy-' ' £(£+1) ^ 

where a dot stands for a derivative with respect to the conformal time 77 and the coupling coefficients are 



(£2_to2)(£2_s2) 



^2 



(10) 

(11) 
(12) 

(13) 



^ In order to avoid confusion with the notations of Ref. [24j . from now on, if not oxpHcitly written, we are going to drop the superscript 
denoting the order of the perturbations. 
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The fluctuation sources are given by 



(0) 



T& 



(0) 



„(o) 



fc*. 



SI 



(1) 



TV 



(1) 



4"^ 



5. 



(2) 



rpio) 



(14) 



with 



pirn) ^ J_ 



(15) 



The modes with m — —\m\ satisfy the same equations with JB|~'™'^ — — i?^^™'^ and all the other quantities unchanged. 
In Eq. (|14p is the baryon velocity perturbation which in the tight coupling regime corresponds to the photon-baryon 
velocity v. 

These equations can be formally integrated, leading to simple expressions in terms of an integral along the line-of- 
sight p3| ■ The temperature fluctuations are given by 



e<"'^(ry,fc) 
2i + i 



dr/e-^5(r^)(V)jf™)(fc(^-^')) 



(16) 



where j^^ are given in Ref. [2^ . For the polarization, we have 



21+1 



Ve 



where the radial functions read 



,(±2) 



(x) 



(x) 



V(^-l)(^ + 2) 

nix)+f;{x) 



(17) 



-Pt'\x) 



x 



-^{e -!)(£ + 2) ^ 



(18) 



(19) 



The differential optical depth f = neaxa sets the collision rate in conformal time, with He the mean free electron 
density and ctt the Thomson cross section and T{rio,ri) = J^" T{rj')dri' the optical depth between rj and the present 

time. The combination 5(77) = re""^ is the visibility function and it expresses the probability that a photon last 
scattered between drj of 77 (with J drjg^rj) — 1) and hence it is sharply peaked at the last scattering epoch. In early 
reionization models, a second peak is also present at more recent times. 



Scalar modes do not contribute to B-polarization, thus 



(0) 



0. We are interested in the contribution to the 



angular power-spectrum for the E and B modes arising from vector (m — 1) and tensor (m — 2) perturbations. 



C, 



(E) _ 



dk 



2 



dk 



m=-2 
2 



= - 1^ l^k 

m=-2 



(2^+1)2 
(2£+l)2 



(20) 



In the following we will focus only on the -B-mode of CMB polarization. This is probably the most interesting 
case when confronting the signal generated by the primordial gravitational waves with that coming from second-order 
perturbations. Both are expected to be very low, but for different reasons: the former because of the low content of 
primordial gravitational waves, the latter because it is a second-order effect. Then, as first pointed out in Ref. p^ . 
there could be the possibility that for a low value of primordial tensor to scalar ratio r, this kind of secondary signal 
could constitute an additional barrier to our ability in detecting the signature of the inflationary gravitational waves. 
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B. B-mode polarization from the tight coupling quadrupole 



The first step is to compute the polarization source term p(™) in Eq. (|17p . In the tight couphng limit we deduce 



from Eq. pip that Gj™' = — 4i?2™VV6i where we have used Eq. Plugging this result back into Eq. yields 



c>(™) 

pirn) ^ ^ (21) 

The B-mode polarization field is then obtained from Eq. ([17]) exploiting the fact that the visibility function is sharply 
peaked around the epoch of decoupling 77,. We thus evaluate P^™'\rii) and p'^p {k{rjQ — -q)) at 77* taking them out of 
the integrals 

= V.)) , (22) 

where we have used A(x, n,?/) = 40(x, n,?/). Notice the different approximations used in deriving Eq. (|22p with 
respect to the procedure of Ref . jl5l | ; in particular we do not need to go beyond the leading order in the tight coupling 
approximation, and this allows an easier evaluation of the time integral (|17p . 

The quantities are related to the quadrupole defined in Eq. ([T|). If in Fourier space for a given k we choose 
63 = k and ei and 62 forming an orthonormal basis, it is easy to verify that 

A20 = -yea.eaj-n!/', (23) 

A2±i = ±5 V3 63. Iflipll n;^- , (24) 

A... = _5y|(fl^(£l^n^ (25) 

where one has to use Eq.([3]) and the explicit expression of the spherical harmonics in the chosen coordinate system. 
Being interested in the B-mode polarization in the following we will make use only of the last two expressions, for the 
multipoles ±1 and ±2, corresponding to the vector and tensor parts of the quadrupole moment. 

The (firs-order) velocities appearing in ([7|) are those of the tightly coupled fluid of photons and baryons, and we 
will use the analytical solutions found in Ref. [l^ with 



7 = -*T17^*kH0)sin(/cc,7y)c, (/C«fce9), (26) 



k' 9 
¥10 
and 

vl^ = -tj^<^'^\0)sin{kcsv)cs (A:»fce9), (27) 

for modes entering the horizon before or after the equality epoch, respectively. Here ^^^"^ (0) is the primordial linear 
gravitational potential on superhorizon scales deep in the radiation dominated epoch, while Cs = (3(1 -I- R„))~^^^ is 
the sound speed of the photon-baryon fluid, with = 3pi,/{4p^) the ratio of the baryon to photon energy density at 
the decoupling epoch. 

Using Eq. ([26]) and ([27]) in Eqs. (p4|) and ([25]) the temperature quadrupole anisotropics can be written as a convo- 
lution 

A^rniKv*) = / |^/™(ki,k-ki)$«(0)a>Wk,(0), (28) 
where the kernel for the vector part, m = ±1, reads 

ki 



- ■ - 20 

/±i(ki,k-ki)k{k^) = T^^c^sin(/ciCs??*)sin(|k-ki|c^?7*)eT''^''i sinOk, 

v6 



2 cos 6'ki ■ 



|k-ki| *^Mk-ki 



,(29) 



a vector normalized to imity, e.g. and that for the tensor part, m = ±2, is given by 



/±2(ki,k-ki) = --?|iVc2sin(fcic,r,,)sin(|k-ki|c,77,)eT2»0kisin0ki^j^^^. (30) 
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To compute Eqs. ((29)) and pO)) one makes use of the orthonormality of the basis = k, ei and 62; (jiki and 
are the azimuthal and polar angles of ki in such a coordinate system. The coefficient N comes from the prefactors 
appearing in the expressions for the linear velocities ([26]) and ([27]) . with N — (9/10)^ if both ki and |k — ki| are less 
that keq, N = (9/2)^ if both the wavenumbers are greater than k^q, while in the mixed cases N = (9/10) x (9/2). 
The power spectrum of the quadrupole anisotropics is defined by 



(A2„(k)A;„(p)) = (2^)3p^,„(fc)<5(3)(k_p)^ 
and, using the general expression ([^5]) it turns out to be 



(2^ 



|/,„(ki,k-ki)|2p^(fci)F*(|k-ki)|) 



(31) 



(32) 



where P^(ki) is the power spectrum of the primordial gravitational potential. The power spectrum of the quadrupole 
anisotropics enters in the quantity we are interested in, namely the contribution to the angular power spectrum of 
the B-mode from tensor and vector perturbations as given by Eq. ([^0]) . Inserting Eq. wc find 



{B){m) ^ 2 f ^1.1.2 \Bi""\k,7Jo)\^ 

7: J (2^+1)2 



647r 



dkk^P^,^{k,r^,) \p^f'\k{7^^-ij,)) 



(33) 



Let us now write down explicitly the power spectra ([32]) using Eqs. ([29|) - (|30|) . For the vector perturbation mode 
we find 



100 r°° 

PA2±iik,V*) = 1/ d(cos6'ki)(l-cos2 6'ki) / dfci iV^fc^ sin2(fcic,?7,) sin^(|k - ki 

3 TT^ 7-1 Jo 



lk-kll 



2 cos 6'ki 



fci 



k-ki| 



P$(fci)P$ ./fc2 + fc2 - 2kki cosOk, 



(34) 



while for the tensor mode 



100 f^^ f° 
PA2±2{k,V*) = — ;^ J ^ d(cos6iki)(l -cos2 6'kJ 



dki 



kt 



sin^(fciCs77*) sin^(|k — kijcgr;*) 



P$(fci)PcE, [Jk^ + kf- 2kki cos 6'ki • 



An efficient way to handle the above integrals over the wavenumber ki is to introduce the variables 



lk-kll 



y 



(35) 



(36) 



fc ' " k ■ 

With such a substitution we arrive at a compact expression for the angular power spectra of the B-mode polarization 



C, 



(S)(±l) 



25 1 
144^ 



dkk^ 



dyy 



dx 



1 



l+y2 _^2^ 2 



2y 



X (j/2 _ x'^y N'' sin'^ (kcsTi^y) sin'' {kcsri^,x)P,s>iky)P<s,{kx) pli^^\k{r]o - 77*)) 



1 2 



(37) 



C, 



(B)(±2) 



25 1 

144^ 



dkk^ 



dyy^ 



11+^1 dx 



i-v\ 



1 - 



1 1 9 ? \ 2 

1 -\- — 
2^ 



X N'' sin"' (kcsTj^y) sin'' {kcsr]^x)P,s,{ky)P^{kx) (3i^^\k{r]o — rj^)) 



(38) 



We remind the reader that N = (9/10)2 if both ki = ky and |k — ki| = kx are less that k^g, N = (9/2)2 jj- i,q^i^ 
the wavenumbers are greater than kgg, while in the mixed cases N — (9/10) x (9/2). The power spectrum of the 
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primordial gravitational potential can be written as P$ = Po*^ '^(fc/^o)"' where Us is the scalar spectral index 
and the amplitude is related to that of the comoving curvature perturbation spectrum, A^(fc) = A^(fco)(fc//co)"°~^, 
by Po* = (4/9)27r2A^(fco). Therefore Eqs. ^ and §^ become 



C, 



(S)(±l) 



25 1 f dk f k 



0*" 



1 - 



144 7r3 y A: V^o 



2^ 



2(ris-l) ^oo 



dy y 



7ls —3 



11-2/1 



sin^ (kcsTj^y) sin'^(fccs?7*x) /J^"''"^'' (^(r^o - ??*)) 



and 



C, 



(S)(±2) 



Pr 



25 1 r dk f k 



0*" 



144 TT^ y fc V fco 



dx X 



ris — 5 



1 - 



1 I 2 2 

1 + y — x 
27 



sin2(fccs?7*2/) sin"'(fccsry*a;) [3''^^' [kirjo - rj^)) 



(±2), 



To perform the integration in fc it is convenient to use the WKB approximated expressions 

4 \ 



1 / 7^2372 



l + ^(£+l)/2 



1 



V.t2 - ^2 



X\J X^ — P 



(39) 



(40) 



(41) 



for X > ^, and for x < £. 

Notice that up to now we have not accounted for the diffusion damping effects in writing these integrals. This can 
be done by replacing the sines with 



sin2(fccs?7*a;) — > siT?{kcsri*x)e ''o , 
sin2(fcc<;77*?/) ~> sin2(fccs?7,y)e ''o . 
where /cd is evaluated at the decoupling epoch and defines the damping scale [27j . 



(42) 
(43) 



IV. DISCUSSION AND CONCLUSIONS 



We show in Figure 1 the results for the vector and tensor contribution to the B-mode polarization multipoles for 
a ACDM model with scalar spectral index n = 0.95 and normalization A^(fco) = 2.3 x 10~^ at fco — 0.002/Mpc 
(the damping scale has been chosen so that fc^^yo — 1600). The tensor modes provide the dominant contribution for 
t < 100, while the vector modes contribution dominates at larger values of £. Notice that the structure of the peaks 
of the second-order effect is richer than the usual (linear) polarization power spectra because the quadrupole signal 
is modulated as a velocity squared giving rise to a sin2(fccs77*) dependence, which appears with two different periods 
in the convolution products 

Fig. [T] reveals that the second-order quadrupole B-mode polarization computed here is larger for £ > 10 than that 
due to the second-order vector and tensor modes, which are produced by the non-linear evolution of primordial scalar 
perturbations, considered in Ref. [l5|. One could also ask what is the analogous signature on the B-mode CMB 
polarization induced by the second-order quadrupole moment of the photons during the reionization epoch. In this 
case the magnitude of the effect has been already computed in Ref. [lg|. A comparison of our results in Fig. [T]to 
those shown in Fig. 9 of Ref. [l^ shows that the contribution of the second-order quadrupole from recombination 
overcomes the analogous signal generated during the reionization epoch by two orders of magnitude in the relevant 
range 100 < £ < 1000. 

It is also interesting to compare this second-order quadrupole contribution to the B-mode polarization with that 
expected from the primordial gravitational waves background generated during inflation. This depends of course on 
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FIG. 1: Angular power-spectrum of B polarization originated at recombination by the second-order quadrupole. The thin solid 
line represents the vector modes contribution and the dashed line the tensor modes. The thick solid line correspond to the 
total contribution. 



the model of inflation considered and its amplitude is parametrized by the ratio of the tensor (h) to scal ar p ower 
spectrum amplitudes r — A^j(fco)/A^(fco) (accordingly to the definition used in the WMAP analysis [H, H^). In 
Figure 2 we show the primordial gravitational wave B polarization multipoles for a model with r = 10~^ (computed 
with the help of the CMBFAST code |30|] ) , along with the total (vector plus tensor) contribution from the second-order 
quadrupole at recombination. We have considered for the ACDM model the cosmological parameters il/^ = 0.732, 
fim = 0.268, Qbh"^ = 0.0218, h = 0.7, t = 0.073 [2§|. The contribution from the primordial gravitational waves scales 
linearly with r, thus we see that for values of r smaller than a few x 10~^ the second-order quadrupole contribution 
becomes dominant for £ > 10, limiting our ability to detect the primordial gravitational-wave background. The peak 
from the primordial gravitational waves at the lowest multipoles, arising from the reionization epoch, could be affected 
by the contribution studied here, only for extremely low values of the scalar to tensor ratio (r < 10~^). 

The main background for the detection of the primordial gravitational waves through their B-mode polarization 
signal is however due to the weak gravitational lensing conversion of the dominant E modes polarization to B modes 
0. Fortunately, this effect can be, at least partially, accounted for through reconstruction of the gravitational 
lensing potential by the correlations of B-polarization between large and small angular scales, which primordial 
gravitational waves do not produce. "Cleaning" of the gravitational lensing signature may achieve factors of 40 
in the power-spectrum, or even larger pT| . In Figure 2 we display the predicted gravitational lensing signal in B 
polarization calculated with CMBFAST for the cosmological model under consideration, reduced by a factor of 40 
[id . ITl| . Clearly, an improvement in the cleaning algorithm would make the second-order quadrupole polarization 
derived here the actual limitation to the detection of primordial gravitational waves through B-polarization of the 
CMB ff r < few xlO-^. 
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FIG. 2: Angular power spectrum of the total (vector plus tensor) B polarization generated at recombination by the second-order 
quadrupole (thick solid line) compared with that induced by primordial gravitational waves (dashed line) with tensor to scalar 
ratio r = 10~^ in a flat ACDM model. The thin solid line is the signal due to gravitational lensing cleaned by a factor 40. 
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